MAURCIO TAKASHI COUTINHO WATANABE - exerccio 5

getwd()
setwd("C:/Users/Maurcio/Documents/aula5")

### Editando alguns parmetros grficos
esaligna <- read.csv(file="esaligna.csv", header=TRUE, sep=",", as.is=TRUE)
head(esaligna)
par(cex=2)
plot(esaligna$ht~esaligna$dap, xlab="DAP (cm)", bty="l", tcl=0.3, ylab="Altura (m)", main="Relao DAP x Altura em rvores de Eucalyptus saligna")




### Dois grficos juntos
par(mfrow=c(1,2))
boxplot(esaligna$dap~esaligna$talhao)
text(0.6,22.5,"A")
# Agora o segundo grfico
media.dap = aggregate(esaligna$dap,by=esaligna[3],FUN=mean) 
media.dap$dp = tapply(esaligna$dap, esaligna$talhao, FUN=sd)
media.dap
plot(media.dap$x~media.dap$talhao, xlim=c(13,33), ylim=c(0,25))
arrows(media.dap$talhao, media.dap$x, media.dap$talhao, c(media.dap$x+media.dap$dp, media.dap$x-media.dap$dp) , angle = 90, code = 2, length = 0.06)
text(14,23,"B")




### Adivinhando o cdigo
# grfico a
codigo <- read.table("exercicio3.csv", header=TRUE, sep=",")
codigo
par(mfrow=c(1,2),bty="l", tcl=0.3)
linear <- lm(codigo$x1~codigo$y1)
plot(codigo$y1~codigo$x1, xlim=c(0.5,2.3), pch=17, ylim=c(0,3), xlab="Log(Patch size)(ha)",ylab="")
abline(linear)
text(2.3,3,"a")
mtext("Euclidean distances", side=2, padj=-3, cex=1.2)
# grfico b
boxplot(codigo$x2~codigo$y2,names=rep(x="", time=6),xlab="",ylab="",ylim=c(0,3), outline=FALSE)
mtext(c("Small", "Medium\nEdge", "Medium\nInferior", "Large\nEdge", "Large\nInferior", "Control"), side=1, line=1, at=c(1,2,3,4,5,6), cex=.7)
text(1,3,"*",cex=1.5)
text(2,3,"*",cex=1.5)
text(3,3,"**",cex=1.5)
text(4,3,"*",cex=1.5)
text(5,3,"***",cex=1.5)
text(6.5,3,"b", cex=1.5)


